Density of states in SF bilayers with arbitrary strength of magnetic scattering 
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We developed the self-consistent method for the calculation of the density of states N(e) in the 
SF bilayers. It based on the quasi-classical Usadel equations and takes into account the suppression 
of superconductivity in the S layer due to the proximity effect with the F metal, as well as existing 
mechanisms of the spin dependent electron scattering. We demonstrate that the increase of the spin 
orbit or spin flip electron scattering rates results in completely different transformations of N(e) at 
the free F layer interface. The developed formalism has been applied for the interpretation of the 
available experimental data. 



It is well known that superconductivity can be induced 
into a non-superconducting material from a superconduc- 
tor due to proximity effect. In the superconductor (S) - 
normal metal (N) bilayers thus induced minigap and the 
shape of the density of states (DOS) at free N metal in- 
terface N(e) depends on the values of the suppression 
parameters of SN interface and relation between the N 
layer thickness and the decay length of the N metal pj- 
[4j. The existence of the minigap has been confirmed 
experimentally in a variety of the proximity SN systems 
(see e.g. HQ and references therein). In the super- 
conductor - ferromagnet (F) bilayers there are additional 
bulk F-layer parameters which influence on N(e). They 
are the exchange field H and the electron spin scatter- 
ing processes. The exchange field tends to align all elec- 
tron spins along the field axis. It splits the minigap and 
density of states for spin up and spin down electrons 
(see H-HIl for the reviews). The experimental study of 
the proximity effect in SF systems [l^- [T^| and Joseph- 
son effect in SFS junctions 0]>E3 reveals that besides 
the exchange field the additional pair-breaking magnetic 
mechanism, namely, spin dependent electron scattering, 
should be taken into account for the data interpretation. 

There are three types of the spin dependent electron 
scattering in the ferromagnet - spin-orbit interaction and 
spin-flip processes which may happen along the exchange 
field direction and in the plane perpendicular to it. Previ- 
ously the influence of the parallel spin-flip and spin-orbit 
scattering mechanisms on N(e) had been considered in 
some limiting cases (rigid boundary conditions at SF in- 
terface, limits of large or small F layer thickness) [Tfjl l20| . 

In this paper we for the first time developed the self- 
consistent method for the calculation of the density of 
states in SF bilayers. It is based on the quasi-classical Us- 
adel equations and takes into account the suppression of 
the superconductivity in the S layer due to the proximity 
effect with the F metal, as well as all three mechanisms of 
the spin dependent electron scattering. We have demon- 



strated that the developed formalism can be applied for 
understanding the N(e) data obtained in superconduc- 
tors with the antiferromagnet ordering. We consider the 
SF bilayer consisting of two dirty metals. They are a 
superconductor of the thickness d s and a thin ferromag- 
net df adjoined at x — 0. All physical quantities depend 
on coordinate x perpendicular to the SF boundary. The 
exchange field is parallel to the SF interface plane. DOS 
can be calculated from the Usadel equations. To pro- 
ceed further it is convenient to use the parametrization 
G(u>, x) = cos 9(u>, x), F(u>, x) = sin 6(u>, x), where G and 
F are normal and anomalous Green functions. The mag- 
netic and spin-orbit scattering mix up the up and down 
spin states which couples the Usadel equations for the 
Green functions with the opposite spin directions. In the 
F layer (x < 0) it gives the system of the two equations 
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and in the S layer (x > 0) the Usadel equations stay 
uncoupled 
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hi(2) = A(x)cos0 s i (2) , (2) 



where 6\ and 9 2 correspond to the Green functions with 
the opposite spin directions, u> — ttT(2ii + 1 ) are the Mat- 
subara frequencies, D s (Df) is the diffusion coefficient in 
S (F) layer, H is the exchange field energy in F layer, 
A (a;) is the superconducting energy gap which is zero in 
F layer. Here we use the self-consistent method to re- 
solve the Usadel equations which takes into account the 
decrease of the energy gap A in the S layer from its bulk 
value along x-axis towards the boundary due to the prox- 
imity effect. The scattering times are labelled here as r z , 
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t x and t so , where T z ( x ) corresponds to the magnetic scat- 
tering parallel (perpendicular) to the quantization axis 
and t so corresponds to the spin-orbit scattering. 

In the S layer the Usadel equations are completed with 
the self-consistency equation 



A(x)\nt + t ^2 



j=0 



2A(x) 



— suit's i — smtf S 2 



= 0, (3) 



where t = T/T c , T c is the bulk superconducting tem- 
perature. Here and further we work with the normalized 
energy parameters A = A/ttT c , uj = lu/ttT c , H = H/nT c , 
and for the length parameters in F layer x = x/£ n , — 
y/Df /2ttT c , and in S layer x = x/£ a , 6 = y / D s /2TrT c . 
The scattering parameter notations are a z = (tzTtTc)" 1 , 
a x = (txkTc)- 1 , a so = (TsoirTc)- 1 . 

The boundary conditions at FS interface have the form 



7s- 



d fl(2) 

dx 
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7 dx 
and at free edges 

98 f 1(2) 

dx 



B =-o= sin(6» s i( 2 ) - Of i(2)), 



U=+o= sin(0 s i(2) - ^/i(2))> 



x=-d f = 0. 



d0sl(2) 

dx 



\x=d,= 0, 



(4) 



(5) 



where 7 = ((T n ^s)/(o's£,n), 0n(s) is the conductivity of the 
F(S) layer, = Rb ^ arL , Rb is the specific resistance of the 
SF interface. 




FIG. 1: The case of the parallel magnetic scattering in the 
S/N bilayer (H/nT c = 0). Spin up energy DOS variation 
in the N layer for several values of the magnetic scattering 
parameter a z and for the fixed 75 = 5, 7 = 0.05, df/£„ = 0.2, 
ds/^s = 10, a x = 0, a so = 0. 

For the arbitrary layers thicknesses, interface param- 
eters, 7, 7s, and magnetic scattering parameters the 



equations Q-JSJ have been solved numerically using the 
self-consistent two step iterative procedure (for ref. see 
[111 13 ) In the first step we calculate the order parame- 
ter coordinate dependence A(x), in the Matsubara tech- 
nique using the self-consistent condition in the S layer. 
Due to the proximity effect A(x) decreases towards the 
SF interface. Then by proceeding to the analytical ex- 
tension in J5J over the energy parameter u> — > —is 
and using A(x) dependence obtained in the previous step 
we find the Green functions by repeating the iterations 
until the convergency is reached. The density of states 
N(e) = iV T (e) + N^e) can be found as 



N ni) (e) = 0.5iY(0)i?ecos6»i (2 



(6) 



where N^n\ is the DOS for the one spin direction and N 
is the total DOS. 




FIG. 2: The case of parallel magnetic scattering. Spin up 
Energy DOS variation in the F layer for several values of the 
magnetic scattering parameter a z and for the fixed H/ttT c = 
0.2, 7s = 5, 7 = 0.05, df/in = 0.2, d s /^ s = 10, a x = 0, 

Ctso = 0. 

The numerically obtained energy dependencies of DOS 
in F layer at the free F boundary are presented on Fig. 
1 1171 At H — (fig. ^) we reproduce the well known 
mini gap existing in SN bilayer 0, Q, Q . The presence 
of the unixial magnetic scattering tends to smooth the 
BCS peaks in the DOS. Figure demonstrates the DOS 
evolution for the spin up electrons for different param- 
eters a z where the full black curve corresponds to the 
usual splitted peaks within the energy gap due to the 
exchange field in the absence of any magnetic scatter- 
ing. By adding the magnetic scattering aligned with the 
exchange field direction one can see the smearing of the 
sharp peaks with the gradual closing of the induced en- 
ergy gap in the F layer. It is interesting to note that the 
symmetry of the spin resolved DOS in respect of Fermi 
energy (e = 0) does not exist in the presence of magnetic 
scattering. 
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FIG. 3: The case of the perpendicular magnetic scattering. 
Spin up energy DOS variation in the F layer for several values 
of the magnetic scattering parameter a x and for the fixed 
H/itT c = 0.2, 7s = 5, 7 = 0.05, d//f„ = 0.2, d s /£s = 10, 
a z = 0, a 3O =0. 
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FIG. 4: The case of the perpendicular magnetic scattering. 
Total DOS energy variation in the F layer for several values 
of the magnetic scattering parameter a x and for the fixed 
H/nT c = 0.2, 7s = 5, 7 = 0.05, d//f„ = 0.2, d s /£, s = 10, 
a z = 0, a so =0. 

Figure |3 demonstrates the influence of the perpendic- 
ular magnetic scattering on the energy DOS variation 
within the energy gap. Total DOS for both spin direc- 
tions for the different values of the perpendicular mag- 
netic scattering parameter a x is plotted in Fig. 0] The 
peaks in DOS are slowly moving towards the zero energy 
that can be explained as the presence of some additional 
splitting field besides the ordinary exchange field in fer- 
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FIG. 5: The case of spin-orbit scattering. Spin up energy 
DOS variation in the F layer for several values of the magnetic 
scattering parameter a ao and for the fixed H/nT c = 0.2, 73 = 
5, 7 = 0.05, df/in = 0.2, d s /£ s = 10, a z = 0, a x =0. Black 
curve corresponds to the SN case. 

romagnet. As in the case of parallel magnetic scattering 
the peaks are smoothed out and the energy gap disap- 
pears. For the small magnetic scattering times t z and t x 
the DOS tends to its bulk value in the ferromagnet. 
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FIG. 6: The case of spin-orbit scattering. Total DOS energy 
variation in the F layer for several values of the magnetic 
scattering parameter a ao and for the fixed H/nT c = 0.2, 73 = 
5, 7 = 0.05, df/$ n = 0.2, d s /£ s = 10, q z = 0, 0^=0. Black 
curve corresponds to the SN case. 

Figures [5] and depict the spin up and total DOS for 
different parameters of spin-orbit scattering, correspond- 
ingly. It can be seen that in contrast to the magnetic 
scattering described above the spin-orbit scattering tends 
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to decrease the effect of the peak splitting within the 
energy gap cased by the ferromagnetic exchange field. 
Black curves in Fig. and correspond to the zero ex- 
change field (SN structure case). The smaller the spin- 
orbit scattering time the closer the curve to the super- 
conductor/normal metal case and two minigap behavior 
degenerate to the one minigap curve as in the SN struc- 
ture. 

It is interesting to mention the peculiarity in DOS de- 
pendence in the presence of the spin-orbit scattering. As 
it was shown in [2l], [22] in the presence of the spin-orbit 
scattering for the parameter r^ 1 = H the solution of the 
Usadel equation changes its characteristic behavior from 
the oscillating one to the damping decay, that should also 
cause the changes in the energy DOS variation. Figure 
demonstrates the appearance of the plateau instead of 
peak in DOS for r^ 1 = H for some parameter 7 when 
it is large enough to diminish the penetration of su- 
perconductivity into the F layer. For the particular set 
of parameters (H,"fB, d s , df) used for calculation of the 
graph in Fig. 0this transformation occurs approximately 
at 7 w 0.5. 
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FIG. 7: The case of spin-orbit scattering. Total DOS en- 
ergy variation in the F layer for r^ 1 = H and for the fixed 
H/irT c = 0.2, 7B = 5, d//£„ = 0.2, d s /^ s = 10, a z = 0, 
a x =0. Black curve corresponds to the SN case with 7 = 0.5. 
Dashed and dotted curves correspond to two different values 
of the 7 parameter. 

Recently, the coexistence of the magnetic and super- 
conducting order in nickel borocarbides was studied in 
several laboratories experimentally. Such compounds as 
ErlSi 'i^B^C and TrciNiiBiC both being the supercon- 
ducting materials demonstrate radically different mag- 
netic properties. Local tunnelling microscopy at low 
temperatures revealed considerable difference in the local 
superconducting density of states behavior. In contrast 
with TmNiiBiC [23 compound where DOS has its usual 
BCS type, ErNi2B-2C 24] measurements show the non 



zero conductance and thereby the non zero DOS within 
the energy gap. 

To find the possible explanation of such a difference 
we propose the following model. We believe that in 
ErN T iiB-iC compound the magnetic order near the sur- 
face is absent even when the antiferromagnetic phase ap- 
pears in the bulk. This may be related with some atomic 
compositional disorder near the surface and modified ex- 
change interaction between magnetic moments near the 
surface. Consequently, to describe the surface properties 
of superconducting ErN 12B2C ^ the model of a thin film 
with the relatively strong magnetic scattering on the top 
of the bulk superconductor without magnetic scattering 
seems to be quite reasonable. 
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FIG. 8: Theoretical fit of the experimental data of 24] (Fig. 1). 
Plot parameters: H = 0, 7 = 1, 7s = 0, d s >• £ S) ct z = a ao = 
0. a) T = 0.15A', T c = 11K, d f = 0.35£„ and magnetic 
scattering parameter a x = 0.95; b) T = Q.8K, T c = W.5K, 
df — 0.6£ n , no magnetic scattering. 

Using the developed algorithm for the SF bilayer, we 
may assume the exchange field H = as in the para- 
magnetic case and 7 = 1, 75 = for the actual ab- 
sence of the boundary. Fig. [SJi demonstrates the cal- 
culated DOS behavior at x = —df in the presence of 
the magnetic scattering which destroys the usual BCS 
behavior. For ErNi2B2C having easy plane magnetic 
anisotropy we take l/r z = and \/t x = l/r y = 1/r. 
FiglSJs corresponds to the case without magnetic scatter- 
ing. It can be seen that both black theoretical curves are 
in a good agreement with the experimental data of [24[ 
(Fig. la and Fig. 2b). The difference between ErNi2B 2 C 
an TmNi2B2C curves may be related with the impor- 
tant difference in their Neel temperatures (6 K and 1.5 
K, respectively). The lower Tjy may lead to the much 
smaller magnetic scattering in TmNi2B2C . 
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In conclusion we demonstrate that the increase of spin 
orbit or spin flip electron scattering rates results in com- 
pletely different transformations of N(e) at free F layer 
interface. The increase of r^ 1 results in the continuous 
suppression of the peaks in the density of states accom- 
panied by the closing of the energy gap. The increase of 
7-jT 1 additionally leads to the shift of the peaks towards 
the zero energy which looks like the action of some ad- 
ditional exchange field in the ferromagnet. Contrary to 
that the increase of t,." 1 does not result in the closing of 
the energy gap and tends to decrease the Zeeman peaks 
splitting. 



All calculations have been done in a self-consistent way 
in the frame of the Usadel equations. The developed for- 
malism has been successfully applied for the interpreta- 
tion of the data obtained in the superconductors with the 
antifcrromagnet ordering. 
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